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partiticns was demonstrated as a precursor for future 
applications to piecewise analysis of oceanic satellite 
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I. ERE ROD UCT ION 


mS M a 
Satellite remote sens 


means of determining the herizontal structure of thé oc¢ans 


@mea globel basis. The objective of this thesis is to con- 


5 Funan dS BS ~ 4 7 ae ~ = : * = ~ he = 
Mmpeere> 2<O =Uhe development of empirical methods £f95r using 


satellite images of optical parameters and s2a surface temp- 


amas ur 


(D 


emeough interpciation and extrapolation of 


meed in situ data. 


(soy tOninte: the Wpeper ocean’s vertical structure, 


relatively lin- 


The objective is being épproached through regional case 


Studies of correlations between cptical para 


ical water mass properties in the upper oc 


regicns cf the world. Mee espe Cas Cally, 


eters and pnys- 


Preemaminaery case study orf the California Current region. 


Men Ultimats goal is to relate statistically the horizontal 


Seeucture ci optical provertiss observed with the Coastal 


Zone Color Scanner (CZCS) to the underlying 
tures of temperature and salinity, és weil 


parameters, for a given region and season. 


vertical struc- 


aS bio-optical 


The study domain enccmpasses the continantal slope and 


emer Off the coast sf California between 


Pee = Arguelic. This aréa was select32i to 


14 


Peene SUS and 


investigate an 





MG etewe eeho S=cnt which is kncwn 5 persist *throughcut 
Bmemtpwel ling season (Traganza,et al., 1979). ie Ce. 2€rn 


PemomecOULRSIn pCrtions of the study domain are typitssa by 
complex addy structure associated with irregular features in 


Sememece nymetrcy, such as off Point Sur. Between Point Sur 


Optical parameters tend to be aligned roughly parailiel to 


the underlying isobaths. 


(D 


An ensemble of data acquired with <=h2 Nimbus 7 CZCS dur- 
=Egeens Summer and fall seascns of 1979, 1980 and 1982 is 
peeeyzeaq In this study. The horizontal structure in o9i0-op- 
Pucal parameters determined from cloud-frea portions of CZCS 
imagery are investigated using a Partitioned Empirical 
Seeenegonal Function (PEOF) decompositicn. Pie "spacial per= 
Mmeeens examined hare consist of four zonal transects cross- 
mmo che Shelf/slope region at differant latitudes. The 
Seec2tic gcais of this analysis are: 

A fommchesace-fizZe tie Merzdional and zonal spatial 
Bercelaticn structures Qi OCsan color parameters 
(specifically optical depth 1/K(490) in méters). 

ae TI¢ compare the Spatial scales and structures of 


optical Vem apa 1 y highlighted by the PEOF 


us) 





Beeompesi ctens, |) end =o relate thsse to the hieterical 
Hesceap-ive Cceano graphy of the study region, and 

Bia To develop preliminary statistics related to the 
feasibility of joining data from different spatial par- 
titicns cn the basis of péertial subsamples, and to thus 
provide an optimal interpolaticn of satellite image 


Gace inte cloud covered areas. 
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iba C@ZA CGR Ooty OCE Tame CENTRAL CALIFORNIA 
To el [SRI eT ee Te 
A. fis SLUDY DOMAIN 
1. Coverade 
a. Area Domain 
The region investigated in this project is 


located between 32 and 4ON, and from the coast of California 
Beeecnore tO approximately 126W, Fig. oe This area was 
selected because it contains water mass structures, includ- 
merrones, Which strongly influence phytoplankton concen- 
*rations, and therefore the optical properties of the océan 
water. Furthermore, an adequate sample orf date was availa- 
Beaees =tCyT this area. 

A subarea of this regicn is labelled Insert A in 
Figure 1 and presented in greater detail in Figure 2 Insert 
Sees founded Ey 34 to 38N, and by 126 to 120w. i 2s tie 
Peete ry study demain of this thesis. Tne background hydrog- 
raphy and dynamics of this region are described in subsec- 


mons 2, 3, and 4 of this chapter. 


ey 





CALIFORNIA COAST BATHYMETRY 


% 


‘Cape Mendocino 





130 W 128W 


126 W 124.W 122 W 120 W 
Pigure 1. Ocean Bathymetry Off the California 
Sas iene: 26 Bazawmetrac Profilin 
stem(SYNBAPS) Data Contoured at 400 a 
erate) 
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bee Tine Domain 


Mime tenmrate:  dOMalim, the 2v¥ailabple CZCS deta 
mmeltuae scenes From summézr through early fall seasons in 


1979, 1980 and 1982. Originally, @ single season ensemble 
of CZCS data (May through September 1980) was sought. How- 
ever, the time span had to be expandsd to three years cto 
assemble a reasonably large sample siz? of cloud-free subd- 
scenes. The sample analyzed contains three scenes in 1979, 
eleven scenés in 1980 and eight scenes in 1982. Detailed 


characteristics of these CZCS images are presented in 


RreweCala fOrnid Coast boundtng tne area cf study 
is characterized by a steep, mountainous coastal range run- 
meme roughly parallel to the ccastiins. iRewececest line 
eeeetchning frem San Francisco to Point Arguello is oriented 
moemagnly northwes= to soutkeast, but is interrupted by Monte- 
rey Bay at 36 45'N and by smallér bays in the vicinity of 
Merec Bay at 35 20'N. Heonmeele jor ELvVersS )dbain intc this 
coastline, although many local rain-generated drainages 


creeks empty here. 


ye 





b. Bathymetry 

WiemeecsCoML Nah. CIZeRtaetion of whe barhyne=ry &s 
meaonly northwest =o southeast (parailel to the coast), Fig. 
2. PAeeEEUDtIONS Of this orientation are evident ir tha 
feeernicy Of the Monterey Canyon, Point Sur and the Sur Carn- 
yon, the Davidson Seamount, the Taney Seamount and the Santa 
Lucia Banks and Escarpment off Point Arguello. A very 
abrupt shelf Ereak is evident all along this section of the 
eae f£Crnia coast. Isobaths tend to diverge south of Monte- 
rey, due to a broadening of the continantal shelf and slope 
Moth distance south of Monterey. 


3. Descriptive Oceancgraphy 


Coastal upwelling is an oceanic phenomenon which 
meena Pronounced impact upen many physical aznd biclogical 
processes. Predominantly southward winds during spring and 
memmer OLL ~he central California coast, yield offshore sur- 
face Ekman transports, which forces compensation water to 


meoes tbom depths of the order cf 200 to 300 m (Smith, 1968). 


is generally confined to the late spring through eariy fall. 


The onset of the seascnal upwelling commences in nore 
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INSERT A 








| (San Francisco 
4200 4 
200 
neon aunt 
4000 . . 
3000 ala 
Se ‘. Point 
15) el onecle ae cael ee 
aes Oeviason 
4400 : Sit yp Seamount 
x 200 
4600 elon | “i “hg er ouetio 
4200 os ; 4 3 
an He : : 
126 W 124 W 122 W 120 W 
meoare 2. Ocean Bathymetry Off the California Coast 
(SYNBAPS Data, Contoured at 200 m Intervals). 
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BSemenernd waters off the California coast an@ ovoregr2ss3s 
Meeenward as «the season unfolds (Yoshida ard Mao, 1357: 
Bemoeo ese: and keid, 1963; Pavicva, 1966; Hickey, 1979). sien 
CZCS data set spans the upwelling season 2nd includes images 
from beyond this season into early winter. 

Upwelling has a marked effect on the sea surface 
temperature, causing it to be much lower than would other- 
wise be normel for the latitude and season (Smith, 1968). 
The relatively lowar temperatures are evident in IR images 
Smmeeie Cegion (Jchnson, 1980; Nestor, 1979). Accompanying 


this decrease in temperature is an incraase in surfac 


(D 


salinity, an upwelling property unique <> the régime off the 
west coast of North América (Smith, 1968). 

The ocsanographic properties of upwelling have 
been documented in many areas of the werld, but nowhere with 
Mmemechoroughness or the work off California and Oregon. 
Ship and satellite observations have allowed us to identify 
Seasons, centers, and the extent of the upwelling event 
meomg ~h2 west coast of Nerth America. hEeada tage, == ai. 
(1979) used combined satellite and shipboard observations to 


! 
Meee nutrient upwelling distributiors off the coast of Cal- 


morn a. Frentzal structures and mesoscale eddies that can 


22 





result from the upwelling phenomenon have been examined with 
ealevance <0 Anti-Submarine Warfare (AS4) Dy »fraganza 
(1979). a ao oe neuased (LR) ahagecy in the @= tection 
and descripticn of upwelling was examined by both Johnson 
(1980) and Nestor (1979). 

PiominlTuaOdieGt tenor nu =rlen= “2¢h waters to the 
nearshore euphotic zone greatly enhances the development of 
fjmenan situ phytcplankton population. This enhancemen= in 
Turn causes the upwelled water mass moving offshore to havé 
Mememence iy different cptical properties than adjoining cff- 
Peete War e€rs. The boundary (frontal region) between the 
upwelled water mass and +he normal surface water mass is 
mms Leacily detectabie and of great interest. 

Hiewehauemenr 1 echmen= Off th] Calafornia coastal 
zone is observed in the régions cf upwelling events. These 
nutrients, which are classified as "bischemicaily newt on 
mem DaSis of nitratée-to-fhosphate ratics which approach 
ior, are brought to the surface from depths up <9 300 n. 
By way of contrast, nutrients also presént in the open ocean 
surface water approach 5:1 (Nestor, 1979). The added 
hitrates are a primary factor in th2 increase in phytoplank- 


*on concentratiors during the upwelling s2ason. COalter 
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Open ocean waters and this dirtference cin be increas? dur- 


ing periods of upwelling (Sverdrup, et al., 1942). The pay- 


iD 


Boplankton concentrations, with ie ae associazed 
@meorophnyli-~like pigments, nave a profound erfect on the 
upwelled radiances measured by the CZCS, as discussed in 
Srapter itt. 

Ani@mnesS aspect ©f upwelling, and its relation +o 
Satellite data, is its effect on regional climates. The reéli- 
atively cold sea surface temperature in upwelling zones 


@mels the air above and thus incr2ases its rélative 


humidity. Neeaetesult, l2Ow stratus and £o0¢ Commonly occur 


{++ 


meer a2 Shallow (marine) layer with warm air alof¢. The 
emeayent occurrence of low stratus and fog, seriously linits 
infrared and visible satellite coverage during the upwelling 
s2ascn. The cool séa water also contributes toa diurnal 
sea breeze by increasing the onshore-offshore pressure gra- 
geen. Onsnore winds bring cool, moist air as far as 50 
Meies inland (Smith, 1968). 
i me CULITents 


The California Current System nay be discussed 


and studied in terms of four large scale currents: 


24 





Bavidcscn Inshore Current, amare sthe SOUgnern Calitorria Cur- 


Bemt (Hickey, 1979). Dicmer: Stu Dhres of * “these Curren<s 
@ecectliy influence the study domain. Mesoscale currents 
associated with seasonal upwelling are also inportant here. 


GQimpice 2 (Or mad CULE Cnt.) fn=  cilitornia Current 
is a broed wind-driven equatorward current which exhibits 


significant seasonal variations proportional to <he change 


in 


mame ne wind field (Brown, 1974). Ort Pomin=, Conmeeo710n -2ae 
Mean annual iccation of the current axis is located 270 ka 
orfshore while the shoreward boundary extends to 200 km cff- 
Or 2. Pacmcreaclcmes Clue Crder 70) xn wide and flows 
south at 10 +0 30 cm per second (Hickey, 1979). 
iecomecalVrcriie “CUSECTt 2S a Continuation of 
ememwest Wind Drirt in the North Pacific and flows southward 
peeeng 2:2 California coast between 48 and 23N. P2=urns 
westward between 20 and 30N where it becomes part of the 
Netenh Equatorial Current. This flow regime comprises thes 
SeWetera ©xtert of the anticyclonic NE Pacific Subtropical 


Gyre, which is centered near the Hawaiian Islands (Sverdrup, 


meee, 1942s Chelton and Davis, 1982). 


ZS 





(Area = tcrmic Undercurrent. Ths Cast Cas 
Geeeoccizcrecnt, ealso, referred te as the California Counter- 
current, is the poleward subsurface flisw over the continen- 


tal slope. Maximum poleward flew occurs during the summer 


fr 


pemgeece lt seasons at devths of 200 to 250 m (Pavlova, 1966 
emeenackey, 1979). The flcw can pe described as a bread 
current with a central jet. Lesa LS ejocm Seructuwe>s thas 
is most often measured and referred to when applying 
specific values *5 the Undercurrent. The broad poleward 
flow has @ ae¢ostrophic ccmponent alongshore near the shelf 
Deak Of approximately 15 cm/sec (Cecddington, 1979). 

The flow appears +9 hava a jet-like struc- 
mee, DOth vertically and horizontally, and to extend to the 
bottom over the slope. The existence of a high speed jet 
Semen ot che order of 20 to 70 km in wiith, was first sug- 
gested by Reid (1962, 1963). Subsequent direct measurements 
of these jets have produced values as high as 40 cn/sec off 
Northern Baja and values of 16 cm/sec off Washington 
(Wooster and Jones, 1970). The depth of the high-speed ccre 
varies seasonally. It rises from depths of 200 to 300 m to 
“he Surface during the late fall @nd wint2r north of Point 


Conception. Here it is referred as the Davidson Inshore 
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@ee-ene DY Hany euthors (Hickey, 19793; Pavlova, 1966; Ingra- 
mama, 1967). evene— Scale Siew uations (Ommene order of 100 
km and 10 days) in the flew appear to be correlated with the 
alongshore component of wind stress (Nelson, 1977). 

The extent and time scale of continuous 
alongshore flow, and the width of the region of northward 
flow below 500 m, are important topics yet to be answered 
about the California Undercurrent (Hickey, 1979). 

(3) Davidson Inshore Current. North cf Point 
Beaception, the poleward surface flow in the nearshore 
ztegicns off z~he West Coast is known as the Davidson Inshore 
Sie rent. It is associated with winter weather circulation 
patterns. As the southward winds weaken and tend toward a 
MOLcThwes-~ward flow, the Davidson Inshore Current becomes 
established (Hickey, 1979). The current flows near the 
eee, USUally within 100 km, well inshore of the California 
Mmerer. and is ccenfined tc the continental shelf and slope. 
Pavicva (1966) reported that north of Point Conception, the 
Davidson Inshore Current reaches its maximum development at 
deeth (200 to 250 m) ain the summer and autumn. In August 


the Davidson Inshore Current is scarcely noticeable at 


*he surface despite active development 2 depth. Maximun 





Surface development is reached from October through April, 


Meee, Le-é€ autumn £0 eariy spring. In December, the cora of 


Maximum velocity emerges at the Surface anc in the lLare 
Spring i+ alnost completely disappears (Reid, 1960; Reid, 
MEmctis, 19508; Pavlova, 1566). Poleward velocities of up ts 
25 cm/sec were recorded (Reid and Swartziose, 196 2) a, eee: 1 
Pemeem Of Central California in January. 

The Davidson Inshore Current and the Cali- 
fornia Undercurrent are often discussed as though they were 
separate currents. Born Clm non ea> eed bslOn- FauactOssal-—*+ ype 
water northward at least as far as Cape Méndicino (Pavlova, 
1966). Also, no suosurface naximum has been found in the 
fiow of the Davidson Inshore Current. These characteristics 
SUpPpcrt a view that the Davidson Current is simply the sur- 
Meee express-on cf the California Undercirrent, rather chan 
a separate current superimrpesed cn it. 

(4) Gener Currents. Th2 presence of eddies 
emeoughouct «he <he California Current System has been decu- 
mented for many years (Bernstein, et al., 1977). The time 
scales for these eddies, as well as the processes resporsi- 
ble fer their generation and subsequent dissipation, is an 


area of increasing study. 
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Between San Francisco and a point about 


tier the distance to Foint Conception, there is a pe 


ry 


mManenc 
@emieesrclockwise eday that produces acrtnward fiow curing 
meeie@on=hs except April (Erown, 1974). A secend eddy just 
north of Point Conception forms during the summer months and 
Makes northward flow continuous frcem Point Conception to San 
Meencisco (Hickey, 1979). Willmott (1983) has shown that 
these features may be preduced by flow separation of the 
Mmm@eeEOriia Current in the vicinity of major coastal capes. 
held, 2t. al. (1963) made direct measurements of 
an eddy (90 km in diameter) Ot Leena HOt cleo upmcons-ecr pata 
Saree crhnia. Hypotheses for eddy formation discussed in 
their paper ate ag Bed GSe 
(1) The process of upwelling and the offshore movement of 
the cclder, more saline waters might degenerate into eddies. 
The lateral shear between the upwelling flow away from the 
@emest and the Califcrnia Current and Undercurrent (dba re— 
mmepic instability) could produce eddy structures. Tempera- 
~ure and salinity differences set up strong baroclinic zones 


along the upwelling boundaries which could result in eddy 


formation. (Sverdrup and Fleming, 1941) 
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OLLeh 


The ore surface flow 


GEN 0ber £Current 


memetan= al north-soutk 


winds, 
Meant occur along the coast. 
(3) 
that the deeper countercurrent 
to the surface layers, at times 
surface current 


Mmmerus:ons of cclder cireulat 


where neither surface 


DEoduce then. 


Malee tionally, Bice ertects - of 


topography, and the associated 


ha) 


Sidered when discussing eddy 


femeson, 1982; Willmott, 1983) 


Cs water Masses 


Descriptions of the 


meomene California Current System 


\ 


(wank, 


VARA 1. On Ln 


A second hypothesis proposed by Reii, 


does not prevail. 
ing waters 


GQOunteccirrents 


+ 


Ou 


Citar WW Slileng, Shot 


195J). 


then seperate countercurrents of different strengths 


Se aiee(loos).,. 2s 
may transfer momentum upward 


when, Or in regions where, a 


This could cause spot 
that form eddies 
nor coastal upwellin 


- 
cd 


Secor img Dy “coastal 
trapped motions must be con- 


Memetiae LOCMmuii Ss. DUE ey, Noy es 


water masses that contribute 


are given in Tibby (1941), 


Sverdrup, ¢t al., (1942) anos ke Ld. “et a). 2 CiS5s)e Four 
major sources are discusse2d by the authors: 
(1) Sukarctic Water Mass - from ths north. 


Centrai Water Mass 


(2) 
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fren the west and northwest. 





ths sewn: 


(3) 


(1941) and Sverdrup, 
USibbarctuc 


These sources were Simplified in Tibby 
naned 


into twe extreme sources 


ae al., (1942) 
Semen Pacific" and "Equatcrial Pacific". 
The percent of each water mass comprising a san- 

= SeOalr. 


eneering figure 3 with a 


ple can ke defined by 
However, the determination of percentage composition by this 
means cannot ke used for water above depths of about 100 on. 
Mme LeSctLictcion is due to vertical mixing in the nearsur- 
face layer related to the effects of wind and local changes 

ane Seats nuetiacs. 


heat and mass fluxes across the 
would be severly 


ct 


due to 
Any mixing along surfaces of constant oa 
masked in these shallower depths by the effects of turbulent 
Also, below 1000 m ths differences in the 


Meecical mixing. 
T-S relaticnships of the two sxtrems water masses are negli- 
‘ONS cle te: 


FOr intermediate depths, as might be 2xpecte 


gible. 
‘ A ® e ® ° . 
percentage of equatorial water decreases in the direction of 
Enarecte -2etically 


Mocthward flow. The Undercurrent is 
Warmer and more saline than the California Current, and Lt 
Gn the oN 26.54 surface. Orr 


has a salinity maximum 
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Monterey and below 3800 


BadetOrial Water and this 


dapth and movement cteward 


Ss 
pee 


SZ 
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Gsagamonowing I-5 Curves Derinin 
augue BeQiiaeoQrial Facific Water, an 
Various Percentages of Eguatorial 
Assuming Mixing Along Surfaces of 
(Beewn, 1974). 
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A. PNTRO DUCT ION 

The physical processes of absorption and scatterinc 
relate the upweliing radiance just beneath the sea surrace 
Mmemmene Constituents of the water (Gordon, 1976). Except for 
coastal waters and waters influenced by river discharge, 
memological constituents play a dcminant role in these pro- 
Mepese= (Smith and Baker, 1978; Jerlov, 1976). Optically, the 
Meee 2mportant biological constituent is phytoplankton, 
Micrescopic plant organisms that photosynthesize and make up 
Mme tirst link of the oceanic food wab (Steele, 197 0) 
Beeorcophyll-a' is the dominant photosynthetic pigment, — 
absorks light Strongly in the blue and red régions of the 
visible spectrum (400 +0 700 nm) (Hovis, ele: VSO). 
Therefor2, as the concentration of phytoplankton increases, 
Semen color of the wacter is shifted toward green hues frem the 
deep blue of its pure state. By measuring upwelled radiance 
(backscaztezred daylignt) in specific spectral bands, we can 
determine the concentrations of phytoplankton pigments in 


Peemocean (Gordon, et al., 1980; Gordon, et al., 1983). 
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[Mis Chavter first déesezibes the CZCS sensor and its 
@eepebiiities, and then the measured signal is discussed. 
Mmfoeicnns that are currnetly applied to this signal to cor- 
rect for atmospheric effects are discussed. Binaliv, tas 
algorithms designed to convert #he corrected radiance values 
memmonycoplanktcn concentrations, C, and irradiance 2ttenua- 


meen co2fficient, kK, are preserted. 


5. mYolLEM DESCRIPTION 
fe che Nambus-7 Coastal Zone Color Scann (S:ze5) 

The cZCS was built by *he Ball Brothers Research 
Mmemepctet2on te NASA's specifications. Pic 8s, rumen. 2S 4 
Spatially imaging multispectral scanner. Six spectral bands 
are precisely coregistered and internally calibrated. The 
Semen Width cf «he CZCS is Slightly mor chan 1600 kn. 
Mmeraecterisctics of its five visible (443, 520, 550, 670, 750 
nm) and one thermal IR (10.5 to 12.5 um) channels are summa- 
rized in Table I. The CZCS has an active scan of 78 degrees 
centered cn nadir and a field of view of 0.0485 degrees, 
Meeeedong a geometric instantaneous field of view of 825 o 
(at nadir) from a spacecraft altitude of 955 km. Ee Gan 


tilt the scan plane 20 degrees from nadir in 2 degree incre- 


fees 2long the satellite track +o mininize the influence of 
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TABLE 1 


Ghamaceerist:c¢cs Of the CZCS 


(Hovis, e+ al., 1980) 


Band Serdnation Measured 
Number Wavelength Gain Radiance Signal/noise 
(nm) (mW/CM5Sr um) 

3 5.41 

1 433 to 453 Z 7.64 1188/1 
1 Bh Ae 
0) 11.46 
3 Bra 0 

2 107 cto. 530 2 De 00 200/1 
1 620 
Q 7.64 
3 2-86 

3} 540 to 560 Z 4.14 1176/1 
1 oe, 
0 Goa 
3 1.34 

4 660 to 680 2 ee 11871 
1 Zase 
0 1, Bisie 

700 to 800 23590 S071 

6 nO, 500 to “12,500 OR 2a k= 


* Noise equivalent temperature differenc2 at 270 K. 
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The designed purpose of the CZCS experiment wes to 


M 
7) 
O 
bh 


provide ¢stinat PRomence=saeface COoncsntrations cf phy- 
toplankton pigments (defined to be chlorophyll-a and its 
associated degradation products, called "phacopiaqments") by 
measuring «he spectral radiance backscattered out of the 
M@eean (Gordon and Clark, 1$81). The radiance scattared out 
of the ocean that reaches the sensor is a very small portion 
of the total radiance received. Consider the physical set- 
sing where solar irradiance POA) at a wavelength Ais inci- 
dent on the top of the atmosphere at a zenith angle 6 and 
azimuth | and the scanner Memo moc ngpto tah radiance L.A) 
ae a nadir angle @ and azimuth angie @. L(A) EQNSIUS=ES Cf 
radiance which has been scattered by the atmospheres and sea 
surface, radiance generated by Fresn2l reflection of the 
direct (unscattered) solar irradiance from the rough ocean 
surface (sun glint), and solar irradiance scattered from 
beneath the sea surface tHE OAD, where t fA) is the diffuse 
transmittance. 

Observations (Gordon, et al., 1983; Gordon, et al., 
1980) produce values cf LOA) in the blue that are ten times 


greater than lives (A) . These effects are principally due to 


oT 





Peacte= ing by the Aix Giewieagh sca tering) aha by micre- 
Peep.c particles suspended in the air (aerosol scattering), 
Moen ct which incrsase the radiance detected at the sensor. 
Presnell reflecticn (fun glint) can be ignored as the tilting 
Seperility or the CZCS minimizes its effect. However the 
scattsring effects, both Rayleigh and aerosol, must be 
removed fren L. (A) to give usable values for the upwelled 


radiance? LA). 


Meee CZCS CEOPHYSICAL ALGORITHMS 


Pactherhes signal descripticn in section Aor this 


Chapter, we can construct the follewing formula 


\ (1) 
LCA) = LOA) + LOA) + tgQOL OA) 


where 
L. = Total radiance 
Pee= Redvance due to Rayleigh scattering 
= = Radiance due to aerosol scattering 
ie eUpMelied radiance from beneath the sea surface 


2 = Titfuse transmittance of the atmosphere 


A = Wavelength 
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+ 


As oreviously nentioned 1b Womens cOtal “seciaxzes 


eee nce CLES. Tre Rayleigh scattering term can be 


Dee tT OA) 


Pao mado! (2!) oh 0(8,)1P (a4) | von 
C2) 


4 TT cos 8 


= The instantaneous extraterrestial solar 


1 


rradiance. 


The Rayleigh optical thickness of the atmosphere. 


The Rayleigh scattering phase function. 


= The scattering angle through which photons are 
Peckscat<-cred frem the atmosphere tc the sensor 
\ 

Veencieo =neeracting With the sea sutfiace. 

Keech Waba | SGCettering angle of these phctons 
which are scattered in the atmosphere toward the sea 
surfece (sky radiance) and then specularly reflected 
BEem ene SUEEaCe into the f£icld of view of the ser- 
SOr (0 (8 .) term) as well as photons which are first 
specularly reflected from the sea surface and then 
scattered by the atmosphere into the field of view 


of the sensor (o9(8) term). 
= The Fresnel reflectance cf the air-sea interface. 


3S, 





ieee =e the fwWe-Wa y OZONE transmittance of the atnoeo: 


O43 ; 

@ = The sensor zantith angle at the observed point on 
the sea surface. 

6 = The solar zenith angle at the observed point on 


the sé€a suriace. 
The aerosol scattering term is found using the 670 


nm channel, where there is only a negligibly small contribu- 


mon by the LS term. CMs LS ESke=rsd =O as the "black 
ocean" assumpticn.) We calculate asrosol radiance at A= 
670 as 

L (670) = L, (670) - L (670) ; (3) 


\ 
The key asSumption in this algorithm is that the ratio of 


aerosol wavelengths 1S constant over a scene, and is given 
as 
LOA) 
e(’,670) = —— (4) 
i a) 7 
a 0) 
Eis calculated using eitkter simultaneous direct radiance 
measurement from a ship, or upwelled radiance values mnod- 


¢lied az a clear water pixel (Gordcn and Clark, 1981). The 


datter method is discussed in section 2 of this chapter. 
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pene nant o seGQudcien (1), ¢. (i) cS ene Gi eTruce 





d 
M@enmsmtitance cf the atmcsphere and sea surface, waich may 
be aprroximated as 
[l ~ 9(8)]} cy) 
Bei A) = —————— exp] - ee 7 Gon mcos oe : (5) 
d 2 03 
m & 


Beeeec=ts heve t2en ereviously defined axcepvt a, which 
Mmemecne index of refracticn of water relative +o air and is 
assumed to be 4/3 for the waveleagths (490 - 700 mm). 

We have now developed the basics for extracting the 


upwelled radiance values, 1,, from the CzcS detected signai, 


The scene constant,e , given in equation (4) is cal- 


Slated as 


Ew@ espa CA) (secs + sec) ] 


Meee) = sal * = - » (6) 
BOA) Saye Oe) + sec) | 
Mere f 2S called the Angstrom coefficient. EaGuar son 9(6) 
can ke rewritten as 
. Fo OA) Bons? 
E(A,A.) = CALA) SSS (7) 
ey) Le ae 
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where Ty (A) = exP| - tT, {A ) (sec 0 + secd. ) and where for 
3 3 . 


mime CZCS ee 670 am. 


(iD 
fu 
(t+ 
aS IE 
M 


Seordon and Clark (1961) develop Goncep. of 
Meer water radiance for atmospheric correction sf cC2zcs 


Mmmagety. The strategy employed in this study was to find an 


area cf the image that cculd be assumed to have a chloro- 


3 
phyll concentration less than 0.25 ng/n. wees Ow cOl 
@entrcation, Li 3t S20m ane 550 8m ate assumed to be 
essentially constant for a given solar elevation. Maen, 


eee nese “clear water values" cf L. (i) at one position, 
{ 

LO) me Celculated using equation (1). et r7o 70) e252 “tcund 

from equaticns (4) using the computed L(A) and L,(670) 


Valus. Finally, cearranging equaticn (7) we find that 


| aa @, aed (dr) 
le (X,A.) J 5 a3 
me A) = 

Ee Coa) Fa hoe . S 





Poe) A) 
0 


Wadues for n at 520 and 550 nm (n(520) and n(550)) are con- 
puted, then averaged to estimate n(443). The Angstrom coef- 


ficient at 443 nm cannot be directly msasured in chis way 
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p2cause2 L442) Peay wocisde ive ~O SVE aL nutes fluccusa- 

Meee Chlorophyll concentrations at low concen*traction. 
miroir eaneaspeece Of this algorithm is that nei- 

+her surface measurements of LQ), nor any properties of 


M@eeee=s=tcsol are required to implement the atmospheric 


aoe entolsophyl Conesrtraticns 
Det ealamcetonmet en oropivyil concentrations 


meccom catios of LA) relates the surface value of © to the 


meetc OL the upwelled Gadianece at two dizferent wavelengtas 
Meee) and Prieur, 1977: Gordon and Clark, 1981). The basis 
Mememeenis is that to a first approximation i is proportional 
memcene Catio of the volume backscattering ccefficient, 

Peeere (A), and the volume absorption coefficient, a(A), of 
m@emwadeet plus its constituents (Gordcn, et al., 1983). The 
Semtributions from the individual constituents can be summed 
memeecrovide a total value for each optical coefficient. 
Memeeover the contributions to B(A), b(A), and a(A) arising 
from phytoplankton and their pigments are assumed to be pro- 
@eeetonal tc chlorophyll concentration C. Taking a ratio of 
L._at+two diiferent wavelengths and applying the assumption 


W 


oe) a B(A)D(A)/a(A), we obtain 
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bCAL) 2 BCA,) BOA{) AQ) 


i Oy) — BONO a Ae 


Meceauss of the non-linearities fever =a ate ne 2 os vyidual 


(D> 


memstatient ccntributions to B{A), D(A), and a(A), we appeal 
moa hewristic cbhservation that 
r 
L., | 1) = )34 Gl Gere meme) ’ (10) 
L Oho) 


meee, che ratio of two upwelled radiances is a function R 
meee Chlororhyll concentration, C, and th2 diffuse attenu- 
meron cOocfficient, k, as well as cther optically important 
constituents of seawater. It was then assumed that 8 is 
Memeered =O C through a log linear medel of the fora 
fog Lo? A. + A; Log R(A1,A9) ; (11) 
which was empirically fit to observations to deternine ccef- 
ficients A. and Aj. Thus, pigment concentrations are conm- 
puted from CZCS data using the equation 
A} 
C=AR (12) 
.e) 
The empirical coefficients oresented by Gordon, 


2t al. (1983) have been adepted by NASA and are: 
Meee t:; C< 1.5 for R (443,550) 


leelizod ooo 


H = 
il 


$4 





Case If: Camis. LOL Perot 550), bus 


Geis. Loe Pogo, os So) 
A eee Vo 7 SOU 
Aye = 1.71 


Meese lil: Cc > 1.5 for Rilo 5550) 
me = 3.3265955 


A, = -2.44 


3 
where C is in ng/m. 


Bemeebaiccuse Attenuation Ccefficient 
A similar de¢eveslopment of the algoritnm for the 
Meee naticna of the diffuse attenuation coefficient, kK), 
memgiven by Atstin (1981). pie meene cnt oroony |? concentre- 
tion algorithm, this algcrithm derives 2 value based on tne 
mr. c Of Le at two wavelengths. k can be defined as 


Raat! aU Os 0 
Ea, Zz) az 


ae (13) 


Memtetszon (13) can be solved for irradiance F(,,z) <0 obtain 


Piezo = F(A,Z)) exp [ -k0X) (eg-23) | : (14) 


Hence 





Jt Kee) 
SS ara 1 FOGzD | ales 
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k(490) = 0.0833 | ——— 40,022 m (16) 


Mestin, 1981). 


Many factors have been accounted for with these algor- 
ithms by either mathematical and empirical models or heuris- 
tic assumptions. The determinaticn of the total radiance 
values in the first four channels of the CZCS allows us to 
apply the corrections to determine upwelled radiance. The 
Bemctituents cf the water which affect its absorption and 
scattering properties are then empirically derived. 

Mie d2stributior of phytoplankten is controlled by many 
local, mesoscale and global factors, including solar radia- 
tion, glcbal weather patterns, anceeecean CazCulation pat 
terns. The mesoscale events OL we Lag or eddy 
Seeeulation can have important regional 2ffects. These fac- 
zO=S are tco numerous and varied te be modellei on a theo- 
mer2cal basis. However, empirical moislling can produces 
relatively accurate and ccnsistent results. 

The measurement of tkese bio-optical parameters from 


Seace elicws us =o remotely determine their relationships to 
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the known (or hypothesized) relationships between inherent 
Mmeaeal properties of the ocean water constituents and the 


Meee involved in their distribution 
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iv. ENZTRICAL QRTHOSONAL EYNCTION ANALYSIS 
he) 


A. INTRODUCTION 

mae, Concept of principal component analysis has been 
Meecenced and utilized in different forms over the past 
eighty years. Freccaigeaesiine tc a data set wes usually 
accomplished using a least squares method. Discances =o 
this line from each point were measured parallel to an arbi- 
trarily set Brie. From the e@arly work of Psarson(1901), 
Mme MSthod was adapted sc that the perpendicular distances 
from each point to the best fit line wers measured. Figure 
4 illustrates this difference and shows that the first 
method is tied ‘*0 a cocrdinate system while the Pearson 
appreach is independent of coordinate systems. This new 
me~hod laid the foundation for the development of principal 
cOompcnent deccmposition teéechnigues. These technigues have 
Since beén utilized in many forms and referred to by similar 
mMames in a number of disciplines. MopinGareons in psychol— 
Semeeey ECkert and Young(1936,1939), although somewhat dif- 
ferent in their development, contain the essanzial elements 
of data analysis and Tincipal component décomposition as 


used in gecphysical disciplines tcday. 
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Axis 2 





megure 4, 


| adie 
eC 1) oS 


Plc. showiag the Difference Betw 
Cem ccanceaetO a Line Parallel) 

Values) and Perpendicular to ch 
(ete fendorkes, et al., 


Seen Zc. On 
to an Axis (d 

Line (p Values). 
Te 0) 
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eewee ete g cal apoliceaticns by Worenz (1956), Coe cr 
meme 196/), and Rinne, et al. (1979) demonstrated the convie- 
Meee Ci representing a large climatclogical fieid with a4 


include: 


an) 
ry) 
4 
(iD 
i) 
tn 


smaller set of values. Thes 


(1) Nen=l2mesr sSta=tiscacGel prediction (Lorenz, 1956), 
(2) Nen-lineaz dynamical prediction (Lorenz, 1956), 


(3) Soo ne Meme c Tela Gomme senvaeion (Rinne, e= al., 


72), and 


(4) Sea level pressure, surface temperature, and ore- 
Cipitation pattern representations (Kutzbach, 1967) 
Other uses of EOF analysis techniqu2s in oceanography 
Mmmenud= the representation of ocsan color spactra (Mueller, 
memroy and of wave spectra (Aranuvachupun and Thorton, 1983). 

Mme Principal difticuities enccuntered in vorincipal con- 
ponent analysis problems relates to the selection of the 
e@eamsangtul’ subtsect cf ccmponents ard *o0 their physical 
Miter pretation. Hemnods@Osesel ection of che principal con- 
POhen=ss are alsc widely varied. Preisendorfer, 2t al. (1981) 
discussed two methods which together involve seventeen Jif- 
Semeemeecesting rules. Empirical selection of a cutoff value 


Meeevyer2zancs or forcing factors can also be utilized. 
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Meeual inerection of the data which leads to a cl2ar cur 


ees an vOptioOn. 


07) 


ferpei= sutjective) choice i 


mabe — =~ 
a 7a SECS 


f}: 


The encermous data volume inherent in satellite 
Megs application of the techniques of principal conporent 
analysis. Frincipal cogponsnt analysis techniques often 


alicw the efficient representation of a large data set by 


ct 


Meetics. few principal ccmponents with 3 negligible loss of 
murormation. The advantage gained is reduction in the nun- 
ber of variables needed tc represent the data. Reducing a 
data set to its principal conponents can also aid in the 
MyeeL pretation cf the data by separating noise from the sig- 
mal. Principal component analysis theory can be applied to 
Meewimi nary explorations within a relatively unstructured 
domain of kncwledge, one in which the fundamental laws gov- 
écning the processes under study are still béing 


defined. (Preisenderfer, at ali., 1981) 


morte f review of the EOP analysis follows to proviie 


background for the later analyses. The reader is referred 
memeerissendozrfer, et al. (1981) for a more complete devel- 
Opment and history. Pie)  LOueeowing) maczix algebra notation 


as adopted throughout this thesis. 
Dome =scomemm@enetes A SCALA. cc ccs cece seo cec kd 


2. A straight line underscore denotes a vector....X 
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GU Neo he UNnGerscote derotss a Bat rlx...s..» 
Pee -naaaqn. line cyverbas denotes a mean values... 


De iim semmorea sclperscrip. "T" denotes a marrix 


transposé. 


Be EOF EQUATIONS 


Following Preisendorfer, et al. (1981) let F' be the 


reg 
va 


Baw (uncentered) data matrix, 


21), 15) se Ci, Dye etae Ue 

zi @2. 1) 2°02) fey ame en 
E' = | 

a ey Cap) 


Meemerr’ (i,j) is the measurement in the i'th time point and 
jeer Spatial positicn. In the present investigation, each 
Meeeec OL Ft will corresfend to an Optical parameter meéeas- 
Meee cy <~<he CZCS at a particular time and spatial position. 
To ccnvert the raw data matrix, Fr, some Centered data 


matrix, F, the temporal means are computed and subtracted 


ecm F*. The temporal mean vector £(x) is calculated as 


Ls 18 
f(x) =n 2 Gees) US 


imencentered data matrix, F, 1s then defined as 


SV 





mG ects’ (le?) — £(2)...£'C1,p) = £(p) 
oem fF 


(imeem — £2) 5.28" (2,p) = Flp) 


apace) © Glos (ne) = £(2)...f'(Cn,p) 


| 
hl 
“~ 
Ze) 
Neer 


Fach elenent f (t,x) ope ey COmsm=e to Orewa caw date Meascurenert 


with the tempcral mean remcved. Whe Centered data Matrix , 


E, man Of written as 


Eee = (20) 


meee F is the matrix containing as rows th 


(D 
ct 
tt 
ev) 
a3 | 
( 
TT 
O 
th 
ip 
O 
tt 


=- - = = 
Smeemean vectcr f= | £01) + F(2) 0 oo Ftd] - 
mee 2Tincipal Dir 


Meommnd 2n2 direction, €1, {in tne s 


‘SO 


ace domain) 
along which the scatter (or variance) Se eeme, Gace. Ser 1s 
Mueemee=-, CCnSsider the projecticn cf the data vectors £ (=) 


mene an arbitrary dirsction $$, 


D ( ie cen? a £ (te, a (21) 
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Mteeong this length and summing over all n observations, 


Be - ~~ - py oy nat aa wn- Vi = be PY om a = ad i ba <5 mm me 
Meta Néeasure Cr tne SCaztser cir the datz along the direc 


Ly 


meon, ©€1, hanely 


2 


n os 
Been) = Dy perce | (22a) 
: 


Mimemrcaghthand side of (22a) can be expanied to yield 


n 
2 
Bece;) = 2 
& = 


a) 
; | = reoetce) [as : (22c) 
ec = 1 


Mimemnext step is to défine the "Scatter Matrix", Ss 


a ee) 2 Cee | (22b) 
l eh ~ ae 


Il 
1 


th 
are. (23) 
Witn elements 
n 
ah 
s(i,j) = DS ere CL) (24) 
Ee el 


Expanding the abcve eguaticn produces for each member of Ss 
Nh 
s eee T 
Meme = 2 (£'(e,i)-F(4))(£'(t,5)-F(4)) (25) 
t 


feeeene Matrix is normalized by dividing by ({(p - 1), then 
when if j, the members of S are covariance? values, and when 
1 = j (the trace of the matrix) «the members are variance 
values (1.€. each element is the variance of ££ ata single 


EMeetal grid tcint). Puc wseaetebelatE x, S, LS Symmetric. 
Ww 
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Mime recOre, i= generally has p non-negative eigenvalues 1. (i 
: 3 
mie o--, Pf) and asscciated sigenvectors Se (Gawler es )s 
pecyided that the rank of $ SseoGuale to p. 
3. Principal Component, Eigenvalue and Eigenvector 
Represenzatcn 
The first principal component of an observation vec- 
met 2S defined to be the linear combination 
if 26 
a; = Sh hla “+ 2 DIED ae) Se ean = - ( 
whose Sample variance 
p Pp 
2 = oe | 
Ss = e,.e,.,.85., = e7,S5 Pad 
a > | > oe eel (27) 
: ios l Fe i 
is 4 maximum for all possitlé vectors ¢, subject to the cecn- 
meeean. chat 
T = 
Intreducing the Lagrange multiplier A,, th] maximum variance 
must satisfy 
3 3 
SE T sb 
S +r, (l-e;,e = e;Se, + A l-e 
ia | (l-e;e)) je, [SS | (l-eye)) 
=e (Sec 47 1) 6) =) 0 ; 
a = (29) 
For nen~trivial solutions,),, must be chosen such that 
(30) 
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A vi 
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8 @1 = Ai &1 Su 
. at 
and since €1¢) = 1 
au 
e S§S e = 4 = s 


fume Sf rst cligenvaluc of S 1S interpretable as the sanm- 
ple variance of S. If we expand this development ‘+o the 


other eigenvalues and eigenvectors of S, we obtain 


E = [3 | ie nero. (. Se) 


the eigenvectcr matrix and 

L = D E | 2) Se aren oe 2 (34) 
mer adiegonal sigenvalue matrix. 
iS » 

SEER (35) 


: aE it 
mmesezms of E micrcOm sera 2nt & 2 = 1 becomes ERI. where 


i Meethe identity natrix. Therefore, if we multiply both 


Sides of equation (35) by Ewe obtain 


S$ =ELE (36) 
NG Nee OG. NG 
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Mmeusing the definition provided for £ and equation (20), 


Mme erincipval component matrix can be defined as 


2 oe (37) 
Thas is the desired principal decomposition of F where 
i = i 

$ a7 (ED EE (38a) 
beeches 
Te lets 8 (38b) 
= = Se 38 
7S % (3 8c) 
= gt EB 2 28 
Gs aa hay (3 3a) 
= | : (28e) 


C. Meer TIONED EOF ANALYSIS 


The EOF analysis method outlined above werks very 
Meme f£o0 a large continuous data sex. However, geophysical 
data sets are rarely continuous. In the case of satellite 
Meee, Cloud cover results in many gaps. Sometimes +hes2 
Meee Can bce bridged by linear interpolaticn, e.g., when they 
See osfall and surrounded by good data. Our enemenis LS hos 
she cas]2 and so a scheme of utilizing non-continuous data is 


hecessary. Here, the purposes of partitioning are: 


Sy 





(1) TO Meximize the sample sige in the presence ofr 
Peod  GOve=, nus aeliowing Stetistical computazriorns for 
eubregions; 

(2) Toma aghlight sratial <ructures of variance fea- 
tures locally, before absorbing them into the modes cf 
the overall domain; and 


(3) To achieve computational convience. 


Briefly, Peele 10nenemepicmemsct¢Ss EOF analysis using 
Small subsets of the overall data set. These subsets are 
Maeo2-.Oned ~c yield continuous data in 3ach subdomain. An 


EOF analysis is completed on each individual subset, and an 
eigenvalus matrix, an Gigenvector matrix and a principal 
eompecnent ma=rix are cbtained. Tee next wseep 1S ace perforn 
an EOF analysis to join the PeameeMale GonpenentsS Ciiac fe sub 
meee his second EOF analysis produces ‘joining functions' 
whach relate twe non-overlapping subsets. 
2. Rules and Methods 

When performing the partiticn of any data set cer- 

fain rules must be cbhserved tc maintain the statistical 


Bemeedi lity cf the computations. Two obvious and basic 


rules are? 
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(2) Peep artl shen sazemsnolldynece bS so smail that the 
Beaeadal Structure is dominated by noise (@.g., a parti- 
tion boundary will not be placed in a major feature, 


er eddy of length scale much less than 


ct 


Such es 4 fron 


The methods involved in the partitioning are subject 
to the abcve principles, tcegether with a general understand- 
mag cf Be physical processes CSE r ngelm che study domain. 
Mepeetrackliines at 35N, 3& 22'N, S35 4O"N end 35 5S3°N, were 
Mumm asd if this initial trial of partitioning (Pig. 5). 
The radiance values and computed cptical parameters along 
Gach tracklinée were plotted versus distances from the coast. 


These plots were aligned to pictcrially reoresent the data 


and its gaps (due to clouds). The partitioning scheme was 
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Mmenme applied £0 try and produce subsets that were 4s cen- 


Meece (COntinuous) as possible over the time domain. Tas 


Meeet data matrix is thus partitioned into >? subdomains 


el. 

Where the surkscript t denctes the <total data set And sub- 

Mmemet p aenctes the partiticns of the data se*+. ach ae 

me data matri PMO Peanea falling WwWLthin gsid parti- 

meme D, and ccntains ali time points for which complete data 
were acquired in that subdcmain. 

3. Equation Development 
iy mEOwmNeecompOosl.2On Glscussed in 8B, is applied +o 


Meee tertition separately, such that for each partition, op, 


the scatzer matrix is given by 
Por! 
and from (36) 


S. = £& L at (41) 
VD VD Pp N 


= Elgenvector matrix of the spatial partition p, 


W 


where 
and 


L = Figenvalue matrix for partition p. 


6 1 





Metco x Of DEincipal componenzs fer each partition is 


given by 

E (42) 
Equations (38) require that 

A , (43) 


mmeso (417) can ke writter as 


S =E ae A gt (44) 
op pe Peep 
Mmeeecing this tc a global scatter macrix, S, 
21 
7 (SYM) 
ca ee 22 (45) 
i 


fee > 1> Lepresents the patzix cf covariances between grid 
Mememes 17 Gdcmains 1 and 2, and so forth. 
Now tke joining process is developed. For any nun- 


Mememcn partitions (two are uséd in this Jievelopment), 


SAG : T i 46 
epee ieee de> Fl Ry Ay Fi Jade) 
and 
ei _ it i 
BO RDN) RCS ea (47) 
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Mummers Cccmbined set the scatter matrix iS giver by 


iE JE 
Fi Fi Eee 
912 = (48) 
Fi Fa #2 £2 
Using equaticns (46) and (47), 
T it sh 
Fi Q AL il aaa Ei 
5912 = 4g 
; 0 ee A ie A ES : | 
Y Fo wl y2 2 2 yY Ge 


Mme) 7CINing Functions, J, are defined as the eigenvectors of 
Mae central matrix given in (49). Finally using (41), (43) 


and (44), (49) becomes 


EB, 0 EL o 
OY av) T av) NG 

a2 = leet tes” 12 T (50) 
AW) av) av) N 


Meee jclaoing functions, oe relate the separate subscenes to 
Memerecr her across an overall study demain. ioc 2nLerpreta— 
Mmem or thes= functions should allow ¢xanination of the var- 
Meereors tha= occur threcughout the domain and Ilccalized 


Meeects On the individual partitions of the domain. 


63 





D. PUL ERPRETATION 


The cocriinate system defined by tne siganvectors gives 
@memdoma:n fer the principal ccmpecnents. In the present 


Slay, as well as most geophysical applications, the princi- 
Mmmeconmponen:s Can be thought of as temporal amplitudes and 
the eigenvectors as their spacial modulators (Preisendorfer, 
eeeeaie, 1981). 

The i-th principal comvonent is that linear combination 
Memone data field which explains the i-th largest portion of 
Mme, total field variance. Essentially the sigenvectors 
define a directicn of variances, Thee oe DEawer bel compe 
mencs give the amplitude cf the varianc2 in the direction of 
the associated eigenvector. 

Once a data sat has Feen reduced to a set of Eigenvec- 
tors and associated amplitudes, the question of signal ver- 
Sus nceise arises. A decision as to which components cf the 
daca field have significance, and which components of the 
data field have no physical meaning must be made. Some sort 
of a selection process must be defined and applied. Bases 
for these selection processes should have their roots in the 


physical processes being studied. 
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A. INTRODUCTION 

Mis Major focus of this thesis was to achieve a first 
Memermtowards the analysis cf the obtained data set. Much of 
the energy in producing these results was directed toward 
Mmenrrocessing of the data to a usable form for EOF analy- 
Sis. Appendices A and B give a detailed discussion of the 
Meocessi 1g <téecnanigues utilized and an accounting of all 
adjustments applied to the data. 

The results cbhtained in this thesis 2ncompass thr2e dis- 
earect areas. The first result e¢merged from the data 
processing and a discovery f the breakdown in the biack 
ccean assumption. The remaining areas are interrelated as 
one deals with the data set prior te tha EOF analysis, whiisa 
Mmiemcthner attempts to relate this to a statistical meaning 


using EOF anélysis methods. 


B. BenxECTIONS FOR NON-ZERO eee) ZN COASTAL WATERS 
Freliminary examinaticn of this data set showed that, 

near the California ccast, the assumption that L (670) = 0 

breaks down (Chaptar III, Section Cc). This finding pre- 


sented a need for an adjustment algorithm. 
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L ($70), Zour ecanety Qreacer than zero, calculated val- 
es ctl L (443) are oe@ten < 0.017 Im Neh Sem (approxi- 
eee 1 )6|6digital)6©6cccunt)6©in CZCS channel 1). Cievse 2S 
unzteasonable €ven in moderately turbid ocean waters. Smith 
and Wilsor (1981) observed that in coastal waters off Cali- 
fornia, where pigments and/for sediment concentration ara 
Mepatively high, it is not unccnmon for the subsurface 
@pwelling radiance L..(670) to be non-zero. They ceveloped 
Mum~meete=atiVe procedure to account for this, which is similar 
*o that develeped independently and used in the present 
PeeccSsing. 
The precedure invelves two major stéaps. 
In the first step, which is invoked when L (44 3) < Ono : 
12 Set I (443) = 0.01, a minimal value for daylight 
Baekecet*e> and sligt<ly less than one digital count in 
SZ7CS channel 1. 
2. Decrease L (670) qicmetersdsc be (070) O= be Consis— 
tent with the new value of L (44 3) (using equations (4) 
and (1)). 


3. Récalculate L (320) and L (550). 


meee Recalculate C,, Co, K(&90Y, K (520). 
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The second step is kased on the assumption that the 
Meme gorithm (equation (12), case III) is robust and insen- 
Sacive +o moderate errors in L (670) . This assumption was 
mmeperced ty sensitivity calculaticns which showed C2 values 
Momvary by less than 30 percent for wile variations in e (A 


Meme rf 6K }lCU WCU lDUC WD COWS 0Cand C1 ~> 2, the correction is pre- 


Sumed unreasonable due to L (670) béing too large (and cor- 


respcndingly, I (670) being t2o small). Next, values of 
K (490) and K(52Z20) are estimated, which are consistent 
Senet Co: 


ae Fstimate a ratic L. (443) /L (550) consistent with 
fe, DY eneene Sos e ce gOr tam (EGUuaelon(12)) with case 
MecoOcificieants) using the C , values. 
mee «6lncrease L (670) and decrease L (670) to be consis- 
tent with the new values of L,(443). 
3. Recalculate a a L (520), and L (350) 
W. Mecalevlare Cy, Cs, K{(4 90), K(520). 
Be) iterate this procedure until C, and C, agree. 
EPaea acguired aboard the R/V Acania during the Optical 
Dynamics Experiment (ODEX) provid2 a tentative basis for 
assessing the validity and performance of the above adjust- 


ment aigorithao. In Figure 6 values of 1/K (490) calculated 
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Memon CZCS Gata, acquired on 16 October 1982, are conpared 
meme D-climinery calculaticns of 1/K(49)) from selscted CDEX 
Staticrs. Macmeranscet Snowe is along 35N  (Sarciticn &). 
mmetrons 24 and 25 were cccupied 1.5 hours before the Nin- 
bus-?7 CZCS observation, and 2 hours after it, respectively. 
Meetion 21 was cccupied S$ hours,and stations 19 and 13 one 
and two days, EeSDeCe LV ely, Gb Ome Oomen> Sceteliite pass. 
The 1/K(490) values at tkese stations were calculated from 
the graphical displays of raw irradiance profiles (at a 
wavelength of 490 nm) presented in the preliminary EF/V 
Serer A ODEX CRUIS&S REPORT Cee, Zaneveld and Smith 
moc2). 

Panel 6a comrares the CZCS and in Situ 1/K (490) values 
before the abceve adjustment was applied, and figure 6b con- 
Bees them after the correction. Agreemsn*+ in both cases is 
eet ient in the transparent waters at stations 21, 24 and 
eee no adjustment for L (670) WesesmcdUaa=2d in this region. 
Ma the inshore portion of the transect, however, agreement 
is obviously poor before tke L (670) aAtgustient, 22a -Dnweh 
imprceved afterwards. This result is pralininary, and sub- 
ject to possible revision by cognizant ODEX investigators 


when their data have been breught to publishable form. 
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Nevertheless, th? EL (670) @davee ner: elastakha so overyh<in- 


megly 2moroves the CZCS estimates of K(990) that its us? in 


Ce 
s 
}+- 
{n 
ct 
r 
D 
OF 
+ be 
tn 
So 
i | 
O 
Jj. 
(iD 
@) 
ct 
pa 
W) 
rh 
ow 
}--! 
}— 
rs 
[LJ }, 
fa 
) 
ct 
} 4. 
rh 
js 
wD 
ju 
jw 
ae } 
ps 
U 
Ww 
tn 
iD 
2 
ct 
\ 
ey) 
ps 


40.0 





2 a. . 
: 2s 7 
he i J F 2 
S | it ie aii! i. 
eee sai: 
Z : \ A | l:> 





: b. a 
2 si saul 
= ’ el i | ! a Nh | a5 
: © y 41} A 2 
: Su ah neil } | 
500 400 300 





DISTANCE (KM) 


Mmegurs 6. Comparison Plots For 1/K (490) Between Track 
Selected ODEX Stations. Panels cy eyigel) \(jo}) 
Hpemelce ve ly are Bertor2 and Arter Ad cmen 
Non-Zero Values of L (670) . 


69 





Meee LATA STRUCTURE 


Figures 7 through 10 shew the cptical depth parameter, 


1/K (490) = Zoo (490) along each track for the available cata 
scenes. (Gordon and McCluney (1975) showed that Zgq (A) es 


the depth over which 90 percent cf La) is backscattered.) 
Mmemerlots are oriented se that the ccast is on the right- 


Mpa side (positive x), while time of the data scene cc2s 


mom earliest to latest in the positive direction along «he 
Beaina-e cf each figure. One scalesogee t(/K (490) as in 
neters. 


miepust Tf and LTII qive background into the oceancgraphy 
Semcene regicn and how that can be related <to ocean optical 
Beare -2rs. The structures depicted in figures 7 through 10 
melt be discussed in terms cf ocean eddy and front visuali- 
zaticns which result from thes2 relationships. Relatively 
high values of 1/K(490) indicates water with lower concentra- 
tions of chlorophyll and sediment. In general these concen- 
traticns may be expected to decrease with distances offshore. 
Abrupt changes in 1/K (490) are usually associated with ocean 


meomcal structure and eddies. 
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= ED ee ee ee ee cee 


The only data available from 1979 (23 Nov) is from 
Mepeet and shows relatively little structure, fies. 1) 
This image was cktained after che end of the upwelling s¢a- 


son and the surface waters were hemogentous to at least 225 
moeotfshore. 

The 1980 data serieS shows mere structure. Begin- 
ning cn 17 May 80, an eddy of approximately 40 km diameter 
Was centered approximately 180 km offshor:2. Sixteen days 
later the entire ‘track shows several eddy-like fsatures 
Mmemnging {Nn size from 4 km to 20 kn. Wasee deys /bee2>) en 6 
meme 50, <he track has lost much of this structure, although 
1/K (490) generally increases in the offshore directicn. 
This trend persisted and strengthened See gqaitly ecisoughs June 
1980. By 1 August 1980, a distinctive pattern had developed 
With rearly unifcrm turbid waters adjacent to the coast, and 
meee orup= (15 km) frental transition to much more transpar- 
2n= waters at a distance approximately 95 km offshore. This 
Pattern is suggestive of thea zonal scal2 of bio-optical 


respense to ccastal upwelling over a single season. 
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The 1982 scenes are ail in the late fall ‘to early 


Oy 


[memees, The transect from 30 Seprepber 1982 shews uniforaly 
Mmeoed water (5 to 10 m cptical depth) Monel etee St uc= 
mire, excep: fora weak, clear-water eddy signature 200 kn 
offshore. Between 30 September and 15 October, more trans- 
Meese (20 to 25 m Optical depth) water intruded from cff- 
Meme tO Within 190 km or the coast, with 42 frontal 
*#>ansizion region of scale appreximately 15 kn. Over the 
ensuing mentn, this clear intrusicn app2ared to evolve into 
meciaecld of less organized, eddy-iike anomalies with scaiés 
meeaagirg from 10 to 50 km. thoeecolo= Mens lone: Character Of 
“has pattern evclution shculd be studied in a future analy- 
Beeemce: chis date set. Such an Ree Which is bsyond tne 
scope c&£ the present thesis, may indicate whether these 
changes are kest interpreted as breakdown of a svatially 
Bees nuous intrusion of offshore water, or as simply due to 
advection Efarsporting “entirely Seeeloic. Waser mass 
Meme ures ints view at this locaticn. 
fee erti tion 2 (Zonal Iransect at 35 40ON) 
rack wo (200 seo) womNOcausd alongs .atecnds 35 


Meme {South of tha first track). 
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Figure 8. The Ortical Derth Parameter, 1/K(u = 
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arene sas 127) data set 25 limited cc one day 


and progresses through late summer (1 August 1980). suey 
Same general features found in the northern track are aiso 
peeaent in track 2d. However, the general appearance of the 
Meenrtel boundary is kroken by additicnal eddy structures. 
Mmeeca:tion, organized océan frontal structures, ww 
apparently related to similar features in Partition 1 data, 
are here displaced approximately 15 km further offshore. 
mime loezesecenss from partition 2 show the same 
frontyeddy develcpment presented in the data from the north- 
Peeeetrlansect (partition 1). Peewee siOs> frontal DO Indas y 
formed between 30 Septenber 1982 and 16 October 1982, and 
evolved +o a less organized pattern of 2idy signatures. 
oo 2 ary2 to! 
Mrondetizeck. 3 (Fig. 9) <the features already dis- 
Sieeecd for the previous partitions occurred. An ad dane so na) i 
feature here is a strong clear-water eddy signature 2aporoxi- 
mately 220 km offshore, evident on the third of Jnne 1980. 


The data suggest an eddy~like intrusion of transparent off- 


shere water, approximately 40 km in diameter. (A feature of 


Us 





Segriar dianeter, eee eee o pam nacely 40 «Kh ciceér 
ashore avoctecars simultanecusly in wvartition 4 (see bealow)). 
Mees teature is short-lived , however, as no évidence of it 
remains *+hree days later. This strange phenomenon is 
another candidate Te cucu re investigation slog a 


2-dimensional analysis. 


In general, frontal boundaries along partition thres 


are farther offshore in both the 1980 and 1982 series than 
apparently related frontal signatures in partitions one and 
TWO. Again, this May suggest the importance of bot«=on 
st¢ering of the mean flow cver diverging isobaths. 
4, Partition 4 (Zonal Transect at 35 90N) 
PH] Seuzhesnmest treck at 35 00 N (Fig. 10) has Sin- 


Mmeepacterins, including the offshore displacement of fron- 


Bee poundaric¢cs. The location of the frontal boundary is at 


Mmeeetarthest cffshore position in this track. 


The transient eddy feature of 3 June 1980, as dis- 
Sassed 1n partition 3, appears here to have a latitudinal 
eeenm: large encugh *o span partitions 3 and &© which are 
Bomugniy 45 km arart. TiicSeomconsas enceweehe che zonal 


scaie of the feature 


(Approximately 40 km). 
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APS ICAL MS 
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Miers al=o iS a longitudinai displacement 
40 km whic Suqges~ssan Gadyeoat Stigntly oblong shape that 
roughly parallels the coast. Questions of whether such fea- 

tures are coherent and continuous from one ‘transect to 
another can ke easily resoived re sueulse ne 2-danensitoral 
analyses. 
It is reasonable *o0 expect a band o£ surface water 
due to 
during tne 


Meet JOW Crtical depth to lie adjacent to the coast, 
productivity 


DCG LegaecaL 
in winten. 


enhanced nearshore 
upwelling season and due tc northward transport of more «ur- 
rent 


Cavidson Inshore 
Ear One t=anspear- 


the 
to be 


by 


masses 
covered by this data set. 


Did water 
Mmeconctcrast, orrshore waters tend 
at least during the seasons 
[Het cn ObVLOUS jont2 cal 
the boundary between what 


EM» 
be anticipated therefore, 
Die -O 0.2 Cal 


It may 

front will persistantly délineate 
and offshore 
is) ot a= 


may be classified as nearshore 
regimes, and that intrusicns of éddy-lik2> surface wate 
che regime to the other will be illuminated by 
Men zeonta, Gradien= 
regime, 


ero ft 
Rodtcacna lly, 
within each 


1, DU eS 
SoeLuctuLe 


Siracal contrast. 

mire] in bio-optical processes will, 

often (but not always) accompany the physical 
ccean frents and eddies and preduce optical 


associated with 


vis) 





Meme st. uCcturs cf Similar scales. The scales of orevisusly 
Mmeerved structures in the Califcrnia Current region may 
thus be expected to be present in the horizontal structure 


of optical depth. It is clear from the above discussion and 


Mmemmiiext Secticn these structures will be discussed frem a 


Breeessticel viewpoint using EOF analysis. 
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Table Il is a listing cz the eigenvalues for tne 
Spatial covariance matrix calculated from che dat2 in each 
Dertition analyzed. The cumulative percentages ef zhe tctal 
Variance ere included. Tre first ten eigenvalues are listed 
Memeen ln each cease, they account for roughly $8 percent of 
*he variance. The €igenvalues are presented graphically in 
figures 11 through 14, which inciude both the eigenvalues 
listed in Table II, and the additional eigenvalues that con- 
tain the nciser-appearing, higher spatial frequencies of 
Seed bat y (and together account for less than 2 percent of 


the total sample variance). We hav2 assumed that 98 percent 


Mmemetne Varlance is an adequate cutoft for calculations. 
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Megure 10. The Optical ae Parameter,1/K (490), Across 
Pidtee seem 4s (35 0.0ON 
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GeemeNejone seaelnies Scoh 25 =GRS tigenvaluss are te 
Mmecetences in tne first value between partitions 1, 2, anda 
See and partition 4. The variance in partition 4 is mere 
evenly distributed over the first thres sigenvalues. ‘Bis 
Meer tat itions kaye reached roughly 90 percent of the vari- 
ance by the sixth eignevalue. The (seers eace 2neiene struc- 
Bare of the eigenvalues suggests that partition four is 
ener affected by additional factors not fourd in the 
Merehe=n partiticns, or that some factors which influence 


the nexthern partitions are absent here. 


S| 


Before the eigenvectors and principal componen<s are 
Mimeeacetreted, how they are combined *o reconstruct a particu- 
lar cbservaticn is explained. Recall mat Seach eigenvector 
Meenes a direction of spatial varibility, and that its 
associated principal comrenents represent the amplitude of 
Meetetions in that directicn at certain time points. In the 
Meeeen. context a "direction" takes ths form of 1/K (490) 
MereetzOns that are coupled at all grid points of the 


domain, SnGeidtrect on’ ~ in this s2nse may be best 


81 





Pertetson 1 Pamet clon 72. . 

Order Eigenvalue Cumulative Eigenvalue Cumulative 

2 2 

(nm )- Percentagé (nm ) Percentage 
1 271.40 4%. 60 Sicpawres ls, 43.70 
2 174.90 68.41 TNS ol soe 
3 LAA 79.47 99.09 81.495 
4 Be i2 Soe. 45.90) 36.94 
5 EAS 5 ShO, Soy ZieaeS Ss S28) 
6 17 369 w2. 0S Z0Ea° 92.83 
44 14.80 ou ,35 zen Vez 94.64 
8 Oye Slash sis Fei! Shomer ad | 
9 Us 97..05 7.90 Sie | 
10 Sree )e: Ws 83 S98 Sh OSA 
Paecitson 3 ele 920, welo)eg" 

Order Eigenvalue Cumulative Bose aN ate Cumulative 

(m )+ Percentage (m ) Percentage 
1 357. 40 Side, ie 273. 10 SE Ais eu 
Z Zo 0 74.65 183.60 Doe > 
5 97.76 Ba. 15 119.80 74.90 
4 a0. OG SEN TES, G4.14 80.63 
5 22.93 SEE SS 4QO.387 §5.93 
6 20. 12 Soe 8, 3 516-96 90.60 
7 13:66 Sis) aye Zee 2 93.77 
8 9.48 Oia 1330 Slo med 
9 8.29 98.0€ Bie2z SG.o 0 
10 Bu 74 98.64 7.00 97 .49 

melustrated either as a GuUave: (£30 one-dimensional 


transects) or centour plets (fcr two-dimensional domains). 
Seven) a temporal mean value for each grid voint of a deonain, 
mmeneigenvector multiplied by the principal component fer a 
Specific time yields a modifier tc the mean signal. 
Poeeoxiple, tO View the Contributions of the first 


five eigenvectors to the observed signal for 3 June 1980 


82 





meaure 11. 
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Eigenvalues for Partition One 
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Beeguite 12. 
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Regure 13. 
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Figure 14. Eigenvalues for Partition Four 
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we muitipiy eack eigenvector By 


The stepwis? reconstructior Se ee observed vec 


atte 3 June 1980 Gata (partition S\ee 2 2On 2S principal 


eigenvectors fe cast rates (eae, 2} 


that the first sigenvector enest (nicer ol ae Sy 


cunts for variations in zonal spatial =e pe hens US i= hale ta 


Bee ot order 200 km ct greater, and with significant 


t 


bitudes only beyon This aspect of tne 


4 80 km of Eshore. 


sle-information contained in eigenvectors aft oe  oebehe 


,sed in later sections. 
1S5b illustrates ihe previously 


mee with @ solid curve, and now <hé Gorted curve 


-ddition cf the second 2igenvecto- as modi- 


dog) by the eecond principal componen- This mode accounts 


; be yariations across thé entire transect with scales Of 


45c shows =e contribution ee ete! 


-der 100 km. Eo ouce 


sird eigenvector and principal ccomponent. This mode has 


tmost no effect O21 +he signal (45 SPSL IE. eh subsequently 4257 


uesed, «his happeas +o be an anomalous Case) -» Piceon ts. 


utions of *ta€ fourth and fifth eigenvectors and principal 


‘ompcnent are ists ate o in Figs. 19q and 19e- Comparing 


fo. 9 with the final A@ctted curve shows «hat “he data 
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ry 


Meera co ttion chree or VIS113 on 25 Jun 1980 and ths recen- 


structed curve aqree with the use of five 2igenvalues, and 


ray 


Mimet a fair tepresentation can be reconstructed using only 


The mean optical depth (1/k(490) = Zoo) transect 
Meeeties £LOr partitions 1 through 4 are illustrated in Fidqs. 
16 through 19 and are repeated for ease of compariscn in 
Figs. cue ~hrough 27. The mean vector in e¢ach track 
represents the tendency of the signal, while the eigenvec- 
tors scaled by the principal compenents give the perturba- 
tions of the mean. Pima leetOUn  Dabule wos, the mean value 
Beeopti:cal depth <tcends tc increase with distance offshore. 
This tendency is expected since the coastal waters should 
Semicain higher concentzaticns of sediment and phytoplankton, 
especially during the upwelling season (Traganza, Gee tar leoee 
1979). There is a general lack of significant eddy-like 
Structure in the mé@an vectors from all four transects 
(although very lew amplitude perturbations of scale five-to- 


ten km are aprarent in the means). 
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meoure 15. 


200 : 100 Land 


Distance Offshore km 


Recorstruction of Optical Depth Transect of 3 
TUN oOMPattition. 3) fren po clmm@so..d Curve, 
Panel a) and Successive Contributions of 
egecivecetors, | to 5 (Dashed Cuzves) in Panel a 
moO Gl ASSpeceively. 
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The eigenvecter discusSicn invoives many intercon- 
parisons of the partiticns. Eacle Sparel tLONs <11S2 Ten 
megenvectors are plotted, (Figs. 16 through 19 and 24 through 
yy. TO Organize the discussion, the first eiaqenvector will 
memecescussed for all four partiticns bafors proceeding ¢ 
d2scuss the second, and se forth. The associated principal 
compcnen<s dae alco tyWistsacea,. (Figs. 20 through 23 and 28 
Paeough 31). 

MeemoerUctGGe sth tho f£LESt Cigenvector of each of 
Mme partitions is characterized by a band of low variability 
Memgpacent tc <he coast, and the structure ofZshcere of that 
band is dcminated by a scale extending from there to the 
Seecncre end cf the domain. Puen Hodge eniarking the onshore 
Timit of significant variation in this mode is progressively 
Mmernet Offshore, proceeding from the south chrough the par- 
meerons, The “node" of partition 1 begins at approximately 
> km offshore, and by partition 4, the "node" is 100 kn 
offshore. There 1s 2 tendency for variance to decrease in 


the amplitude cf the first ¢€igenvector as the offshore 


boundary is approached. This may be an artifact of the 
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Smter bOUndery and should be investigated further cver a 
Meege: dOrain to better estimate the dominant scais. Parti- 
#ion 4, which has the larcest spatial extent, shows more and 
larger offshore structure to beycnd 180 km. 

The associated principal components, Vinci Od ane ¥ 
“he e€igenvectcrs before they are afpliei to the mean, show 
the timé variations. Mewessrcne four tracks, whe, figs. 
eigenvectcrs/principal components vary in phase with each 
other. In all four transects there is a large difference 
beztween the first principal component of the only 1979 image 
(early winter) and those from 1980 (early spring). This 
marked difference is certainly a manirestation of seasonal 
Meeeeaticns in the California Current system (Pavicova, 1966; 
muey, 1979). Most cf the first principal components vari- 
ability in all cases is okserved in the 1980 series (upwell- 
2ng seascn), and the record contains relatively littis 
Meeeaadility in the 1982 series (Davidson Current season). 
Meerererncy cf the variations differs from partition-to-parti- 
Memon with nce apparent pattern. The first eigenvector and 
principal component appear to haves their foundations in the 
ofishere sé€ascnal variaticn and large scal2 eddy structure 


that cccurs during the upwelling season. In the first mods, 
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Mm/emersnore zone intlinenced by upweiling tends to remain 
merc thrcuchout the year, whereas ths dominant variations 
memortical depth occur offshore of the upwelling zon:. 

The shapes of seccnd esigenvectors from the four par- 
Mimerors are Similar, but very in an oscillatory fasnion fron 
Mmemnorenéern cartition to the southern partition. The first 
Partition (northern) shows negative values beyond approxi-~ 
mately 180 km offshore, and then small amplituds positive 


values from there to the coast. Patti con two -and three 


Mmgeect a Mirror lmage pattern to that of the first parti- 


aren . waGeacaOn LOUS Shews much the sane pattern as parti- 
e2on one. The phase relation in the principal components 
mows 20 pattern between partitions one and two, Dc, eae 


Peete s fCr partitions *+hree and four both Suggest phase 
Meversal from the first partition. Phas Aegative-to-posi- 
tive-back-to-negative pattern cf behavior weakly suggests a 
Mere-iike meridional cscillatory structure, with an cffshore 
meee (180 km cffshore) in the vicinity of partitions two and 
Sees. Resoluticn of this meridional characteristic feature 
Will require a 2-dimensional analysis. The distance of sep- 


aration of the partitions suggests a wavelength of the order 


Seeeizo km. Again, the majority cf the variability occurs 
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meomse 16. Mean and Eigenvectors 1 +0 5 for Partition One. 
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meee 1/. Mean and Eigenvectors 1 +0 5 for Partition Two. 
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Peogure 18. ueee SueEGenmveetons 1 £0 5 E£Or Partition 
Three. 
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Mean and Eigenvectors 1 +05 for Partition 
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PRINCIPAL COMPONENTS 1 - S 
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megane 20. Principal Components 1 to 5 for Partition One. 
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meagre 21. Principal Comronents 1 to 5 for Partition Two. 
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PRINCIPAL COMPONENTS I - S 
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Mere 22. Principal Comrcnents 1 to 5 for Partition Thres¢. 
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PRINCIPRL. COMPONENTS 1 - S 





Time Point 


Mmegmee 23. Principal Comgcnents 1 to 5 for Partition Four. 
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furing the 19680 series. MOBS) N=arSshOr]e Structures are 


meee bilt+y, as compared tc the first E€ligenvector. 

The third eigenvector has a similar behavicr for the 
mest “WO partiticns. The perturbations are of roughly the 
Same spatial scale (45 km) and appear to bé in phase. Hcw- 
ever, she third and fourth pazrtitions show an opposite 
behavior in the far offshere région (beyond 180 km). Numer- 
ous smaller scale features (of the order of 10 km or less) 
are apparent in this eigenvector. In general, higher spa- 
fial frequencies become increasingly important in higher 
ordér eigenvectors. ire eps ncatal “components Show an 
increase in the variability of the 1982 series with «he wide 
Tange of variability still present in the 1980 series. This 
eigenvector shows the largest nearshore amplitudes of all 
the eigenvectcrs, which suggests it may be closély linked to 
meme mearshor=s structure cf upwelling. There is little sug- 
gestion Sect temporal pelacion evident in the third princi- 
Beeeccmponents of the four tracks. 

The fourth eigenvector shows an increase in fre- 
quency (decrease in wavelength) of the represented variabil- 


ity scales. Feétures range in size from 18 - 45 kn, with 
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numerous smaller scalé perturbations. nie éaveragj2 
Wavelength and tange of variability is approxinstaiy «he 
meme fOD all four partiticns. Picmore nei S41, COMPORSh<cs fOr 
partitionsone, three and four qualitatively suggest coher- 
ancy. However, @igenvector behavicr is opposit2 in parti- 
Mion three and similar in partition four when compared <«o 
Maetition one. A prominent feature approximately 180 ka 
offshcre in the fourth eigenvector of partition one, seems 
memoe shifted cutward to 210 km offshore in partition four. 
This time/space relationship between the structures of par- 
tition one and four again suggests a meridional oscillation 
worthy of future investigation through a 2-dimensional 
analysis. 

The fifth eigenvector shows features Of SCake se hat 
range from 5-42 kn. The much smaller features (less then 5 
km) are not dealt with as they are essentially part of the 
background ncise expected in any natural system. Little can 
be said akout the correspcndence of the four partitions with 
just a visual inspection. However, nearshore structure 
appears in ‘the eigenvectcrs with more variability than in 
previcus eigenvectors of the 1982 series. The 1980 series 


Still demonstrates the largest overall variability, and «he 
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and «he 1980 scénes is still evident. i) sey) (2) eal lees iE am folewnelel 
Meres Of partition one and four show excellent agreasment in 
amplitud2 and phase for the first five time points, but as 
Structural variatility decreases with ‘tine, so dces corre- 
spondence. The quantitative correspondsnce of these varia- 
memes s beyond the sccpe cf this thesis, but it should bs 
investigated in future analyses of this data ser. 

The sixth through the tenth eigenvectors are charac- 
Memmeeeead DY Variation cf such high frequency that little can 
be said of the relationships between the partitions. Scales 


Seeestruccural reatures in these eigenvectors range from 1 to 


iD 


35 kn, With neo suggesticn of a temporal relation between 
Peet tions. Thea peice components show that the 1980 
data again dcminates the variance, but at these higher fre- 
G@aencies the increased contributicn cf the 1982 data to the 
total variance is very apfarent. Eecause Of =nis disorgan- 


Meeqg structure, detailed interpretation is not attempted for 


Seegemvectcrs cf crder greater than 5. 
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meoguse 24. Mean and Figenvecters 6 tc 10 for Partition One. 
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meeagure 26. wee” admmegenveceor =e 6 £e 10° ror Partition 
Three. 
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PRINCIPAL COMPONENTS 6 - 10 


2.8 


10 


Time Point 


Megs 289. Principal Comrcnents 6 to 10 for Partition One. 
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mmore 29. Principal Compenents 6 to 10 for Partition Two. 
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mure 30. eee cevad Compenents 6 te 10 for Partition 
hree. 
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Sena JCTe2bg Of Imic Pereiticrs 
Pie emeamelnoonye thee more Of che joining function 
results is presented. Partition one and two were ezxanined 
and analyzed és per the jeining frunction development given 


an Chapt 


analysis on partitions one and 


cde 
Ne 


paring 
analyses 


manct*cn 


computed 


more, using ten degrees of freedon, 


and 
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two combined, and tren con- 


+o these obtained fron 


7OLnC a. ieewecerOund merhat th] jJOlning 
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e 


Mesen ing the Srincipal components of the jolming func- 


joining precess. 

Bae e2genvalues OFteained by thus joining =he eligen- 
Mectors and principal comegcnents frem partitions 1 and 2 are 
Given in Table III, together with fractions of total sample 
Paciance. At fac? value, the first ten eigenvalues account 
meg e.595% of total sample variance. Recall however, <zhat 
the input data were represented in truncated form, using 
May the first ten principal components from each of the 
Moeecitions. This original appreximation retained only 
Bieso9e Of the total variance computed ‘from the originai 
data, hence, it is nécessary to adjust the apparent trunca- 
meron cit the joined result accordingly. The results of this 
adjustment are given in the third column of Table III and 
Show that assuming that only the first ten eigenvectors are 
Begmeracant actualiy leads to a ‘“*runcation «=o 96.50% of the 
total sample variance. While 3.5% precision is an accepta- 


Mpe level cf approximation for mest problems in g¢sophysical 
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Meeee Interpretation, pe eh tees O8 §SUCeessSiVs =SLncaticn 
Mee be given carerul attention when applying the parti- 
tioned methcd to EOF analysi 

This example indicates strongly that <he pazrtition- 
mag approach to EOF analysis can provids computationally 
acceptable results when applied tc satellite image data. 
The exampis also amphasizes that proper care must be given 
Mumm Olling Successive truncation in the partition join- 
ing process. ooo re co ri uire Projects to investigats 
@uestions suck as joining partitions on the basis of par- 
emally intersecting Samples to provide optimal functions 


meme ncetpoOleting satellite data inte cloudy reqions, and + 


O 


Mipeetpre.-4ticCn of partiticn joining functions to illuminate 
Spatial ccrrelations between lecally important structures 
(@.g., tcpographically generated mesoscal2 eddies) and the 
dominant structure of the cverall domain (e.g., that associ- 
ated with the evolution cf the synoptic scale upwelling 
frent over the continental slope and shelf over the course 
femeene upwelling seascen). It is questions of this kind «hat 
address the ultimate utility of partitioned EOF analysis. 
The present effort is linited tc preliminary werk to estab- 


Mmeomesoundaticns of feasibility and procedural constraints. 
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Zonal transects of optical depth (t/k(49G) mm) measured 


with the Nimbus-7 CZCS have been analysed +o investigate 


Meee so. . Chl Structure sver the ccntinental shelf and Siopea 
Meeecentral California. Samples or cloud free data were 


Mmmected and processed fer latitudes 35-53N, 35-40N, 35-22N 
and 35-00N. The data were observed in 1979, 1989, and 1982 
during the months May through November. The zonal structure 
in these samples was analysed using EOF's computed sepa- 
Bately for each section. Me¢ridional variance structure was 


m@eueeeysed Only qualitatively through inspection orf similari- 


ct 


Mmee in tsatures contained in EOF's of he different «rarn- 
sects and in the temporal sequences of associated principal 
compcenencts. Bena, Picmeetpucacaonai  f£seSibiiiry of 
memvyang partiticned EOF an@élysis methods to this «ype of 
data was investigated by joining the EOF's of the two north- 
eacnmost <cransects to form estimates of the EOF's of the con- 
bined spatial domain. 


The first eigenvectors for four zonal «ransects of opti- 


cal depth 1/(k490) ¢ach contained dominant scales of order 
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Meoekm OL g=eater, and acccunted for bewween 35 and 54 per- 


Memeeby 4 band cf low variability in optical dep+~h in tke 
Meshore region influenced by upwelling and the Davidson 
Mronere Current. Pies eaves mecOnraned within 45 km of the 
Meameeceat 35 53' N, and menotonicéelly broadens to approx:i- 
feeely 100 kn a* 35N Latitude. This béhavior is possibly 
celated to tne broadening of the continantal shelf and slcpe 
Mich Longshore distance scuth of Mcntersy. jel Bho lLalbbence ff a) feck 


showed that the topographic beta effect plays a fundamental 


role in the dynamics assceciated with mesoscale crdéer 100 
kn) MemgshotTesvariations in tcpegraphy by affecting the 
Besength ct the longshore flow. Alseye <=he inflwence of 


topography can produce barectropic flow beyond its immediate 
mecinity. For mesoscale variations in coastline geometry, 
memcoastal currents and the patterns of vertical motion 
mto LOllow the coastline, but net wich unirtorm strength. 
Coastal current widths tend to be narrower than the scale of 
coastline variability. In these terms, the meridional vari- 
ation in scales present in the first EOF's are qualitatively 
comsisten+t with the longshore variations in bathymetry of 


the study domain. 
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Gominant scale cf order 120 kn, ANG Mee Nearly uritcra 
em@mtbaitude from the coast te a node approximately 150 km cri- 
Mmere 20 partitions 1, 2, and 4. The sscond éE€igenvector for 
Sititaon 3 (35 20N) is ancmalous in that 12 is dominated by 
MeezOnal Wwavetorm with nodes spaced at approximately 80 kn, 
Semecughiy halt the deminant scale of its counterparts. The 
reascn for this behavior should be investigated. 

The chird eigenvectors are dcminated by scal2s rangin 
from approximately 60 to 100 km (between nodes). The shapes 
Mpaeescales vary more strongly frcem partition-to-partition 
than was the case with the first two eigenvectors. 

Across each transect, zonal features with wavelengtnas 
100 km and greater appear. Ths sugg2stion Of anh oscillatory 
behavior in the meridional direction needs to be studied 
mast her. Resclving such a feature requires a more derailed 
study involving a 2-dimensional analysis of the study 
domain. 

The large eddy field associated with the shoreward 
boundary of the California Current was observed in the data 
set. The scales of this eddy field were of the same 


magnitude as the spatial scales employei in the partiticns. 
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bourdery within 


Mey ficiG anc cutting away somes of che feactur3s. The sea- 
sonel development of &@ Synoptl SCatc do veileng ~s20nt off 
the Califcrnia ccast is strongly suggested in the data and 
jts eigenvectors. The sgaaéller eddies 2ssociated with this 
Baectern ranged from 5 to 100 km in scais. 

The convergence of the eigenvalues to roughly 98 percent 
of the variance after the tenth value was of particular 
interest. Dicmewceoenec CUC efOL brail  f£Sur partitions and 
although this is not an overwhelming reduction in the 


degqrees of freedem of the 
Satellite images, and 


pheric variables, provide 


initial systen, 


cther 


massive data sets. 


teers Suga. Can. . 
EX Ses OL Oceanic and a*nces— 


Large amounts 


of computer time must often be expended for processing these 


data sets 
interpretations are, 


volume of data. 


mathematically representing satellite data 


pact and 


EOF's illuminates, and facilitates analysis of, 


at even relatively primitive levels. 
morever, 


EOF analysi 


easily manipulated forn. 


Analyses and 
made difficult by the sheer 
provides a viable method for 
fields in a con- 


Data transformed using 


Space scales associated with a given variable over the 
domain; the present study has exercised this attribute of 


hac 


4 : “a -_ x ‘eo ~ om Seiad en —_- . a ‘ane ae -_ = = 5 Pred o ~~ +78 we So 
M~~m@eoarpa. CCNCONEnt LCepresentaticn Cr Satsiiite mages sre 


j 


eo 
ivt, 


vides an efficien form for analysing the response of spa- 
Meeweece cr uctuse in, for example, optical depth to forcing by 
Wind stress and currents, Aon Ge -NoOUugn, a. -o1O0-c sti cal 
Memeets =his is a logical avenue fer future research to build 
on the present results. 

Considéring purely computational aspects of EOFs, ths 
well-known symmetry of eigenvector solutions in the time and 
Space domain can be used te great advantag2? in the analysis 
of satellite image data. The number cf spatial grid points 
in even the single trackline partitions of the present 
Scudy yield large, but computationally tractable, scatter 
matrices. The larger arrays asscciated with 2-dinensicnal 
area partitions, each with several hundred grid points, will 
clearly exceed sizes admitting direct computation of spatial 
EOFs. The linear algebra and scalings involved in using the 
smaller time domain SCcamimensgtatatk fOr Computation of space 
domain EBOF's is reviewed in Appendix C. 

The partitioned method of EOF analysis illuminates cor- 
relations retween variability in spatially separate sub-re- 


q2ons. The present results demonstrate the comoutational 


Zo 


Mecical depth date. There 1s #svery reason to believe thax 
the method may be equally weil applied t29 other CZCS pérense- 
Be-s and tc infrared imagery of SST. BHimenocewoak, 2n  chis 
acea should ain to determine whether the joining functions 
ienkaing EOF's from separate domains ars sufficiently site- 
MmeonaLy tc provide a basis for optimally interpolating sat- 
ellite image cata of these types over cloud-covered areas of 
@ particular day's image. Other applications to be explorad 
[mmelude determination of the extent to which correlations 
between 3-dimensional in situ data and 2-dimensicnal satel- 
lite data in small sub-regions may be extended to other 


parts of the larger domain covered by satellite data alone. 
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APPENDIX A 


See h eee eee ROG SS ING METHODS 


Ars INTRODOCTICN 

Lata processing was divided into three major levels. 
Mevel-1 processing includes all steps required to ‘take the 
original data tape to a Level-I tare. Level-II processing 
includes all steps between a Level-I tape and a Level-It 
tape. Lével-III processing includes the stenos involved to 
meee the Level-Ii output to a usable forn. The Col lowisg 
sections briefly describe the steps involved in the three 
levels of processing. 

Computer hardware utilized. was that resident at the 
Naval Postgraduate School, Montersy, California. The main 
frame computer used was the IBM 3033AP while the mini-con- 
puter used was the Apple<II. Computer software referred to 
mmeenes SeCtich is either a system utility resident to the 
IBM system or a locally generated program. Documentation of 
the lecaily generated programs can be obtained fron: 

Dr. Jz Le. Mueller (Code 68My) 
Department of Oceanography 
Naval Postgraduate School 


Monterey, Califcrnia 93943 
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programs cn other systems are cautioned to review tne dcecu- 
mentation carefully prior to attempting to *«ransfer the 


Sortware. 


Semeee eve L—-i FEROCESSING 


Figure 3Z is a scktematic diagram illustrating the 


in 


processing foo eroOr elbow eoueammana, should be referred to 
emeQuonrout this discussion. The master tape (raw satellite 


Mmeeame was obtained from the Scripps Institution of Oceancg- 


rh 


saphy, San Diego, California. Table IV gives a summary o 
*he master tapes utilized in this study. The data were in 
Mmemcesm Of aistandard magnetic tape in a binazy format with 
BeeeurbitsS per inch (BPI). The tapes were originally created 
using a Hewlett Packard (HE) - 3000 which has a characteris- 
@eeen-gh crder, low order bit arrangement opposita to the 
IBM systen. PuerelOLre,  Leftors using this raw deta in the 
MemeesUs3AP, it had to undergo a byte swap routine. ian 
byte swap was accomplished when the unforma«téd working bac- 
kup tape was made using local pregram VISBKV. After the 
unformatted backup tape is made @ variable blocked spanned 
(VBS) format tape is preduced using ths system utility 


IEBGENER. 
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Figure 32. Level-I Processing Schematic Diagran 
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Sateilite 
Tape Source 
Pees gnaticn 
mr50 17 lempus 7 (eZCs 
mes 032 Ninbus 7 (CZCS) 
YIS040 Nimbus 7 (CZCS) 
VIsSog4 Nimbus 7 (CZCS) 
mesS095 Neanbas 7 (CZCS) 
VIS097 Nimbus 7 (CZCS) 
VIS104 imbus 7 (CZCS) 
vES105 Nimbus 7 (CZCS) 
VIS106 Nimbus 7 (CZCS) 
VIS117 Nimbus 7 (CZCS) 
VISi11 Nimbus 7 (CZCS) 
yes i i2 Nimbus 7 (CZCS) 
W2S113 Nimbus 7 (CZCS) 
VIS126 Nimbus 7 (CZCS) 
AROO000 Nimbus 7 (CZCS) 
AR2642 Nimbus 7 (CZCS) 
AR2668 Nimbus 7 (CZCS) 
AR2685 Nimbus 7 (CZCS) 
AR 2686 Nimbus 7 (CZCS) 
AR2691 Nimbus 7 (CZCS) 
AR2693 Nimbus 7 (CZCS) 
AR2704 Nimbus 7 (CZCS) 
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*apes serve as the werking tapes for the remainder of she 
Level-I processing, and the master tape is archived. 

Using the VBS formatted tape, a Véersatec plotter grays- 
cale is procuced using local pregram ZIPSIO. This vrogranm 
alse unpacks ‘the Event file (hereafter referred to as the 
E-file) and writes it tc a stcrage disk. The grayscale 
depicts the satellite pass in picture form for hand analysis 
of landmarks. line numbérs and vixel numbers are taken off 


the grayscale for clear, cloud-free landmarks. These valuss 


iw 
tf 
p) 
ct 
iD 
ra) 


are entered into lecal progran F630 21)2) 1G S& aclo mee = tal 
merc PRINT'. This 1S a matrix of radiance value centered on 
she individual landmarks line and pixel number. These 
PICPRINTS are then contoured by hand (uSing a ctzhreshold 
value of 18 ccunts for land or clouds) to determine an exact 
time and pixel number fer the landmark. The landmark's 
latitude, longitude, line number ard pixel number with addi- 


tional housekeeping data are eantered into local progran 


ct 


CZCSNAV on the Apple II. This program is interactive and 
prompts for necessary inputs. Additionally this progran 


adjusts roli, pitch and yaw to reduce the root mean squared 
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Mee ance error in the navigaticn freblen. The msan ins 


YWaiue chtained for all the adjusted, utilized data was 
meprckamately 1.09 n.ni.. i=in SeOaue. Gf ~tThis Ste2o 


generates a set of navigation parameters that are used to 
Meterats 42 navigation matrix. This step is accenplished 
memng local program CZCSNAV2 to generate the navigation 
matrix and NAVDUMP to write the navigation matrix (here- 
after referred to as the G-file) to a temporary formatted 
tape. 

mie E=—-file is copied from its temporary disk storage to 
the temporary tape as the G-file. Additionally the Data 
file (hereafter referred tc as the L-fils) is first reversed 
memes CCT ICM—<O-toOD Orientaticn to a top-to-bottom orien- 
Set2cn using local program BACKWARD. This program also puis 
mae L-file tc the previously mentioned temporary storage 
sape. Finally, these files on the temporary tape are copied 
to @ Léevel-I tape using the system utility IEBGENER. The 
only difference between tke temporary tape and the final 
Memeo -t tape is that the L-file is copied into an unfermat- 


ted file which will aid in the speed of further processing. 


ea 


oe me si-2k 2£sOC=SSi NG 

Peaugee 52 aS a sChémacic diagram illustratizg the 
orocessing steps for Level-Iil. The Level-I taps generated 
by the steps Giscussed in the previous section is the input 
merece fOr this processing. Cmivetne L-s2le is affected by 
mae Level-II processing as the E-file and G-fil=2 are covnied 
mecesgnt tO the Level-II tape using the system utility 
MerGleNeER. The L-file is used to generate output for calcu- 
Mitrng the proper vaiues cf the Angstrom coefficient for 
each scene. Pes seden]e —useng Local program CZPARMS2 and 
an assumed value for the Angstrom cce2fficient. Chapwes Ltt 
Meecren C.2. discusses tke importance and method of finding 


a 


These values. Next, the ccmputed Angstrom coefficients wit 


J 


meme L-Lilé are srerun thircugh CZPARMS2 to regensrate the 
L-file. This regeneration involves taking «he raw counts of 
Seem channel and applying the bio-optic algorithms discussed 
mueenaplier Tit to produce values for chiozophyll and K. A> 
this point the adjustment algorithm discussed in Chapter V 
has not been applied. Analyses cf the initial Level-II cut- 
Pueeecrecipiz~ated the need for the corrective algorithna, 


Which was then applied during Level-III processing. 
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Figure 33. Level-II Procéssing Schematic Diagran 
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D. Poko eit PRCCESS ING 
[ete ee Uoetaeee sta LevVesiciii processing Stsps. 
The Level-III processing basically takes the data ontainsd 


in Level-II and marries it to the navigation aatrix gener- 


mereodgurang Level~-i processing. USe2Gn ene Our dasiquat=< 
meecksc (Fig. Seen e@ "IGOascal Starting points ft>om each 


Bedok were entered into lcecal prcegram TEDDUMP to provide the 
Mmm~gactiton bleck of the track origin. The G-fil¢ contains 
data (in latitude and lcngitude values) every sixteenth 
pixel and sixteenth Line. Once the origin block is estab- 
Mmomeg the exact line and fixel was interpolated using local 
BeoOgnamn FINDPIX on the Apple II. Womens) Startang  pcint 
local program DATAY was entered to generate every 1km along 
each track an associated line and pixel number which was 
shen converted into the appropriate daté values. rons Cue — 
Biemwes Written to storage fer later processing. It was 
here that the adjustment algorithm was applied, producing 


the final version of the data in a navigated form. 
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Pigure 34. Level-III Processing Schematic Diagram 
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peewee kW 5 


DATA CONDITIONING 


Ic apply the Level-III data te the analytic techniques 
Certain conditioning steps were necessary prior 9 begin- 
a9 . Much of the conditioning applied to the data was 
dependent on the data itself as to its completeness and 
DehavVior. This discussion focusses on the steps necessary 
Preto: tO using the EOF analysis techniques. 

Figure 35 depicts the steps involved in this discussion 
and shculd be referred to as a guide. Pincmeciemtey Cl— 11 T 
meee £O> ~26CtEEe)6four.)6tracks and twenty-two scenes were 
extracted using local pregram PACKJOB. Twenty-two files 
ememecontdainec the data for the four tracks for @ach partic- 
ular scene. These data were plotted using Ilecal progran 
PARPLOT and the DESSELA “Weadities resident on the 
mee SO33AP. The format of the plot was chosen *0 give an 
MmmeeecatiOn Of Either good data or bad data with ac struc- 
ture. Mee p! Se Wes wsed to Gecide an the partitioning 
scheme. Four partitions_were selected and their details ars 


listed in Table V. 
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TABLE 


if 
Partition dimensions 


Track Bee Leer on dene Grad Max Grid Ntime Nspace 
No. NO. 
1 1 1&0 410 17 Jogo 
2 z ies 425 17 Zo 
3 2) 200 467 16 268 
ub us 200 485 A 286 
megucres 36 through 39 show the plots generated by par- 
plot and the partitioning given in Table IV. ane Objective 


of the partitioning was tc find the most complete data cover 
time and space possible given ten samples. 

ig@ea crom chesie four partitions was then entered into 
(iereak PLrOgGram PARDAT. This program applied most of the con- 
diticning to the data set. Only the K(490) data was wti- 
Mezeds from this point on althcugh this program could be 
easily altered to focus on another optical parameter. The 
data wer2 search2d to find good points and bad points and a 
control arrangement for later use wes made. The raw data 
were scaled and inverted to preducs 1/K (490) values in 
meters. Next the data were averaged by every fourth peint 
to smooth out noise features. At this point data strings 
with gaps existed for each applicable scene. Next a linear 
Meeempolecting routine was applied to obtain continuous data 


a= Gach =ime point. Fanally, the data for all scenes and 
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MayokS Werte comoined and written to disk in a ¢ingié Gasca 


Mmmeaeel ct oOnzng the data into four partitions, rejecting 
incomplete scenes, Scaling the K (490) values, averaging th2 
data by every four values, and applying linear interpolation 
momtill in tne remaining gaps. 

The conditiored data was then plotted using Ilscal cro- 
gram PARPLT and th2 DISSFIA system utilities. Ene plors 
generated are ficures 7 threugh 10 and were used in Chapter 
¥V Section B to discuss the data and its relationship to the 
regicnal oceanography. 

[ne tinal steps of the data conditioning involved appli- 
GCaticns of the EOF analysis techniques. Local program PEOF! 
produced eigenvalues, eigenvectors, and principle components 
Moepeeeach partition and plotted the output. Pigures 11 
mmeougn 14 and 16 through 31 are the plots produced. The 
eigenvalues, eigenvectors and principle components were all 
Meeeeewen tO disk for later Ned. Biemlocs berpmogran JEOF2 was 
G2signed to produce the jcining function that related par<i- 


meomecne tc partition two. 
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Partitiion Two Date 
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Meee 38. Partitioning Scheme for Track Three (35 22 N) 
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The desired EOF's are those in the space demair, which 
Por satellize data, is dimensioned much lacsger than che time 
dimension. ipeiswpossiole 2.OC Significantly expedit> compu- 
Pareeens by computing eigenvalues, eligenvectozs end principal 
@empenen-S using the smaller covariance matrix cf the time 
Memain, and to then scale these results to oktain the eigen- 
Meesers and principal comrcnents in the space domaizx. The 
algebraic basis fer this approach is raviewsd in this 
appendix. 

Mensa der the followine convention for dimension nota- 
tions 


meacS -« «© « « 2 = lrevesN 


PeIeCOEd=sT) s « «© « K = leeee,K 
where K< min(M,N). As kefore the raw data matrix is given 


by 


ao 2 i (55) 
~ mn = 


of M rows ry N cclumns. The sample mean in the space demain 


is given by 
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Seeaimensicn 1 X N. 


Dnemecen-ered date matrix 
domain is given by 


From Chapter IV she sample space ccvarianca matrix (of 


and the sample tim? covariance natrix (centered 


of size M X M) is now defined as 


(59) 


ic e, are an associated eigenvalues and eigenvector 
respectively cf S, then 
WMA 


— 


= (60) 
B Sy Se Ue 
since 
i # 0 and e, # 0 (61) 
Then 
F ey oie ners (62) 


and therefore 





—— = re 2 
eee ek kw Ek (63) 
(noting the ccnmutative property of hy and & Nise Aste aS 


an ¢igenvalue of T with associated eigenvector (in the time 
N 


domain) 


where a ak ioe tic mm ne Didneipal “Component (in the space 
domain) at time point m. 

We now consider the eicenvalues A and eigenvectors E 
Ok r, normalized such that 


; | 
) wy ont Soe (65) 


m= 1 


with time dcmain principal components 
oe F ex cee Lome cy (66) 


Mmmemerthe ‘significance order K is selected either to retain 
Sumeeesbicrary fraction of the total sample variance, eps 
according to one of the more objective statistical selection 


rules di 


in 


Gucsed Ey Preisendorfer, e«t al. (1981). 
The principal components in the time and space domain 


have the property that 
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Recall that F was centered on € in the space doiean , we a.d 
a = t 
Meee cpcotn S and T were thererore norméiized by tne factor 
ays te 
(1/1). 
Ey analogy tc equation (63) 
Ih a ak Bay 
ss = A. F & =a are 
mre Fe &, ee k 2k? Gee 
where ox 2s now seen to be an eigenvector of S, pues Of 
length 
= rv ~ 
meecmer then kEeing an ort hcnormal eigenvector of length 1. 


The two are thus related by 





1 


fetrneci pal comronents in tke time 
| ag 


mal eigenvectors in the space domain), 


and by similar arguments 


ai) a: ee ek 


(Orthonormal eigenvectors in the time domain 


Cipal components in the space domain). 
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(70) 


dome tieescale tO Orthonor-— 


(71) 


scale to prin- 
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